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ABSTRACT 

Coalescing massive black hole binaries are produced by the mergers of galaxies. The final stages of 
the black hole coalescence produce strong gravitational radiation that can be detected by the space- 
borne LISA. In cases where the black hole merger takes place in the presence of gas and magnetic 
fields, various types of electromagnetic signals may also be produced. Modeling such electromagnetic 
counterparts of the final merger requires evolving the behavior of both gas and fields in the strong-field 
regions around the black holes. We have taken a step towards solving this problem by mapping the 
flow of pressureless matter in the dynamic, 3-D general relativistic spacetime around the merging 
black holes. We find qualitative differences in collision and outfiow speeds, including a signature of 
the merger when the net angular momentum of the matter is low, between the results from single and 
binary black holes, and between nonrotating and rotating holes in binaries. If future magnetohydro- 
dynamic results confirm these differences, it may allow assessment of the properties of the binaries as 
well as yielding an identifiable electromagnetic counterpart to the attendant gravitational wave signal. 
Subject headings: black hole physics — galaxies: nuclei — gravitational waves 
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L INTRODUCTION 

Electromagnetic signatures of the coalescence of two 
supermassive black holes would, in concert with the 
detection of gravitational waves from the event, allow 
a high-precision determination of both the redshift 
and the luminosity distance to the merger. This 
would provide a precise probe of cosmology, limited 
mainly by uncertainties in the magnification, and hence 
the lum inosity distance , due to weak gravitational 
lensing ([Hughes &: Hold l2003f ). As a result, there 
has been substantial recent interest in mechanisms 
that could produce electromag netic emission from a 
surrounding accretion disk (lArmita ge fc Nata,raianl 
2002; Milosavlicvi c fc PhinnevI [2005: Dotti et al.l l2006t 
Kocsis et a l. 2006; Phinney 2007 : IBode fc PhinnevI 

20071: [Shields fc Bon ning 200 8: Lipoai et all l20'08F 

Schnittman fc KrolikI 120081: iKocsis fc Loebl 120081: 



Kocsis et al.' '2008': 'Haiman et al.' '20091 'O'Neil l et al 
2Q09; Haiman ct al. 2009; Phinney 2009; Chang_eLal 
20091 : iMegevand et al. 2009). However, there have been 



relatively few investigations of the electromagnetic sig- 
natures that could emerge from the dynamic spacetime 
near the bina ry during its last few orbits (for a recent 
exception, see iPalenzuela et al1l2009f ). 

There are reasons to believe that if significant gas ex- 
ists around the binary in its last few orbits, then intense 
but short-lived electromagnetic signals will be produced. 
For example, since the binary orbital speed can be half 
the speed of light or more, even in the Newtonian case 
the slingshot mechanism could be particularly effective 
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in ejecting matter or producing collisions that have a 
high Lorentz factor. In addition, it has been proposed 
that magnetic fields in the vicinity of the orbit could be 
wound up and amplified to such a degree that they would 
produce a strong Poynting outfiow (R. Blandford, per- 
sonal communication). Exploration of such effects would 
require simulations that involve magnetohydrodynamics 
(MHD) in a dynamical spacetime. 

Here we take a step towards such simulations by trac- 
ing the paths of pressureless test particles (i.e., parti- 
cles of non-zero rest mass in geodesic motion) in the dy- 
namical spacetime of coalescing black holes. Our goal 
is to determine whether the outfiows and collisions are 
sufficiently high-speed, compared with the equivalents 
around a single black hole, that unique signatures of the 
binary motion seem plausible. This will pave the way 
for the extensive development required to do a full anal- 
ysis of MHD in a dynamical spacetime. In § 2 we dis- 
cuss our simulation setup, which uses a code that has 
been used to model successfully the coalescence of black 
holes of different masses, spin magnitudes, and orienta- 
tions. We choose units in which c — \ and G = 1; with 
this, the total mass M is related to time and distance 
by M ~ 5 X 1Q-^{M/Mq) s - 1.5(M/Mq) km. In § 3 
we present and discuss our results, in which we exam- 
ine first the coalescence of equal-mass nonspinning black 
holes and then show that equal-mass black holes with 
aligned spins produce significantly more extreme effects. 

2. SIMULATION SETUP 

In the following simulations the spacetime was 
computed by solving Einstein's field equations us- 
ing finite differencing methods of num erical relativity 
([Imbiriba et all l2004t [Baker et al.l 120081 ). On the ini- 
tial slice, t he black holes are represented by "punctures" 
((Brandt fc Briigmann 1997) and the constraint equa- 
tions are solved using a multigrid solver ([Brown fc Lowd 
[200l . The black holes are then evolved forward in time 
by integrating the time-dependent Einstein equations 



2 



VAN METER ET AL. 



given by the BSSN formulation (jShibata &: Nakamural 
I1995I: iBaumgarte fc ShapTrol [l998h . We use coordinate 
conditions that a llow the black ho le s to move freely 
across the grid (iBaker et alj l2006t ICampanelli et al.1 
120061 Ivan Meter et al."2006'rMeanwhile the test matter 
is modeled by non-interacting point-particles. To evolve 
the test matter, we integrate the geodesic equation for 
each particle in this spacetime. The metric at the lo- 
cation of each particle is obtained by interpolation from 
the computational grid. 

The geodesic code was tested for a variety of trajecto- 
ries around a single black hole, both spinning and non- 
spinning. For coordinate conditions in which the con- 
version to Boycr-Lindquist coordinates was known, the 
particle worldline was compared with independent cal- 
culations and found to be at least 4th-order convergent 
in the grid-spacing, as were the constants of motion, en- 
ergy {E), angular momentum (L), and Carter's constant 
(Q). For the coordinate conditions closer to those of our 
typical binary black hole runs, the constant of motion L, 
readily calculable in any axisymmetric coordinate sys- 
tem, was still found to be at least 4th-order convergent. 

Initial data for the spacetime background used here 
was chosen to represent two black holes, each of mass 
m — M/2, in nearly circular inspirals. In our simula- 
tions the black holes were either spinless, or were each 
given a spin of a/m = 0.8, aligned with the orbital angu- 
lar momentum. The initial linear momenta of the black 
holes were given by the post-Newtonian approximation, 
and the initial separation of the black holes was chosen 
so as to give at least five orbits before merger. As a con- 
trol case, we have also run equivalent simulations with a 
single isolated nonrotating black hole of mass M. 

In these spacetimes, approximately 75,000 geodesic 
particles were initially distributed uniformly throughout 
a solid annulus of inner radius 8M, outer radius 25M, 
and vertical full thickness lOM. We excluded particles 
within the inner radius to avoid transient signatures from 
particles initially near the horizons. We have chosen such 
a geometrically thick disk because these are the disks 
that potentially have high enough inward radial speeds 
to keep up with the shrinking binary, as opposed to be- 
ing stranded at large radii as the binary coalesces (see 
iMilosavljevic fc Phinnevll2005f ). 

We explored two initial velocity configurations for the 
particles. In the first configuration ("orbital"), the ini- 
tial velocities are randomly distributed around a tan- 
gential velocity Vc that would give a circular orbit in a 
Schwarzschild spacetime of mass M, resulting in a scale 
height of 5M. In the second configuration ("isotropic"), 
we consider an extreme case where the particles only have 
random velocities, with each component sampled from a 
Gaussian distribution of standard deviation Vc/V^- 

We output the positions and proper 4- velocities of the 
particles every ~ IM and ~ 0.5M in the nonrotating 
and rotating cases, respectively. Particles within a hori- 
zon are discarded. We use the relative velocities between 
two "colliding" particles as an energy estimate of the ma- 
terial within the disk. We define a particle collision as 
two particles traveling within Tc — O.IM of each other, 
assuming that they travel on straight lines between suc- 
cessive output times. We have found that the collision 
energies are not sensitive to in the range of O.OIA/ to 
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Fig. 1 . — Maximum Lorentz factors of collisions between particles 
as function of time before merger for initial velocities that include 
an orbital velocity plus a random isotropic component. The three 
panels are for the following cases: (top) isolated nonrotating black 
hole, (middle) two spinless black holes, merging in approximately 
five orbits, (bottom) two black holes each with spin parameters of 
a/m = 0.8, merging in approximately six orbits. In all cases both 
black holes have mass m. The data have been aligned with respect 
to the time axis such that the common apparent horizon is first 
detected at t = 0. M is the total system mass. 

IM. For every particle, we search the nearest 8 neigh- 
bors for any particle that has a closest approach less than 
Tc, which we mark as a collision. 

With this list of collisions, we can now obtain the 
Lorentz factor of each colliding particle, in locally 
Minkowskian coordinates, as follows. Given two colliding 
particles A and B, first we compute the scalar product 
of the proper 4- velocity of particle A with the proper 
4-velocity Ug of particle B. Note that regardless of the 
original coordinates in which the product is computed, 
it exactly equals the negative of the Lorentz factor of 
particle A in the rest-frame of particle B (as it would be 
calculated in locally Minkowskian coordinates). That is. 
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and g'^^ = 77^1/, where the primes re- 
fer to the inertial frame of particle B. We then trans- 
form to the center-of-momentum frame by solving for the 
Lorentz factor ^com characterizing the Lorentz transfor- 
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3. RESULTS AND DISCUSSION 

Our main results are displayed in Figures [1] through |4l 
In Figures [T] and [21 with initially orbital and isotropic 
particle velocities respectively, we show the maximum 
Lorentz factor of a collision between any two particles as 
a function of time, for several different setups including 
one with an isolated Schwarzschild black hole for com- 
parison. Because collisions do not affect particle tra- 
jectories in our calculations (i.e., the geodesic paths are 
unaltered), we expect these to be upper bounds on the 
true Lorentz factors of collisions given that in the MHD 
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Fig. 2. — Same as Figure[T]but for particles whose initial velocity 
is purely random and isotropic, with no net orbital component. 
From this figure we note that (a) in both binary cases there is a 
clear maximum in Lorentz factor just prior to merger, whereas (as 
expected) the isolated case has a steady Lorentz factor, and (b) in 
the spinning case, the collisions throughout and near merger are 
more violent than in the nonrotating runs. 
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Fig. 3. — The fraction of particles with positive radial velocities 
in a shell with surfaces at 25M and 30M with a corresponding 
Lorentz factor greater than (top to bottom) 1.2, 1.4, 1.5, comparing 
the cases with an isolated spinless black hole (solid black) , a binary 
with spinless equal-mass black holes (dotted red) , and a binary with 
equal-mass black holes with spins of a/m = 0.8 (dashed green). On 
the scale plotted, one particle gives a fractional value of 1/75,000 
or 1.3 X 10-5. 

situation realized by nature there will be dissipation of 
relative speeds through collisions. 

In Figures [1] and [H the top panel shows the control 
case, in which the particles orbit around a solitary black 
hole. The other two panels involve nonspinning (middle 
panel) and spinning (bottom panel) binaries. The merger 
in these cases, defined by the detection of a common 
apparent horizon, takes place at t = 0. 

Our control case shows as expected that after a 
short initialization phase the maximum Lorentz factor 
is roughly constant. Here the highest-speed collisions 
are between particles that have both fallen close to the 
horizon and thus have acquired speeds considerably in 
excess of their initial orbital speeds. Given that the par- 
ticles have a "thermal" velocity distribution in which the 
random speed is comparable to the orbital speed, a rea- 
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Fig. 4. — Same as Figure [3] but for the case of particles having 
only initial random velocities and Lorentz factors greater than (top 
to bottom) 1.2, 1.5, 2, 3. Also shown are the cumulative sums of 
these fractions. 



sonable fraction of particles is captured by the black hole. 

In the binary runs the varying accelerations caused by 
the binary motion are expected to alter significantly the 
energy and angular momentum of particles within 2 
binary orbital radii, just as in the Newtonian slingshot 
analogue. More rapid orbits therefore are expected to 
lead to higher-speed collisions and (as is evident in Fig- 
ure [3]) faster outflows, with the effect maximized shortly 
before merger because this is when the orbital speeds are 
greatest. 

Examining Figure[2]we find that these expectations are 
largely met for the isotropic particle velocities. There is 
indeed a strong peak in the maximum collision Lorentz 
factor just before merger. This peak is highest when the 
holes are rotating, as is reasonable because for the aligned 
spin configuration the holes can get closer and thus orbit 
faster before they plunge together. The Lorentz peaks 
are in all cases considerably higher than in the isolated 
hole runs, which might lead us to cautious optimism 
that this phase might be observably different for a bi- 
nary merger than for an isolated hole. However, in the 
case of particles with initial average velocities that are 
approximately orbital, the central region of the disk re- 
mains evacuated near the time of merger, as the particles 
have found quasi-stable orbits further out. In this case 
there is no discernible merger signature, as seen in Fig- 
ure [TJ We conjecture that the inclusion of viscosity would 
push orbiting matter closer towards the binary, resulting 
in energetics between the above two extremes, and recov- 
ering a merger signature. If this is confirmed by future 
detailed MHD calculations, it may provide another way 
to identify an electromagnetic counterpart to the grav- 
itational waves from a binary supermassive black hole 
merger. 

Other features of these figures are, we believe, arti- 
facts of the limitations of our simulation. For example, 
even with isotropic particle velocities, in both the spin- 
ning and the nonspinning binary runs we see an overall 
decrease in the maximum Lorentz factor of collisions dur- 
ing inspiral. We suspect that this is because our setup 
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includes no angular momentum transport between par- 
ticles, hence as the simulation progresses we run out of 
particles that would naturally plunge towards the holes 
and produce very high speed collisions. 

In Figure[3]and[4]wc explore the outflow Lorentz factor, 
measured in a shell of thickness 5M with its inner surface 
equal to initial disk outer radius of 25M. Particularly 
in the case of isotropic particle velocities, because we 
have a significant random speed in our thick disks, the 
isolated hole has some fraction of particles that move 
rapidly outward, but we see that the fraction of such 
particles is (as expected) very small compared to those 
of the binary runs where particles are flung out via a 
slingshot effect. 

The main meaningful difference, in Figure|4l is between 
the spinning and nonspinning binary simulations. At a 
given time before merger, due to spin-orbit coupling in 
the spinning case, the holes are orbiting faster than in 
the nonspinning case. As a result, the slingshot effect is 
stronger at all times from merger, and near merger itself 
the speeds are higher yet. It is also possible that frame- 
dragging due to spin plays a more direct role, by increas- 
ing the ejection speeds in a manner akin to an eggbeater. 
At the highest speeds seen, there is no parallel for the iso- 
lated black hole run. Depending greatly on the nature of 
the interactions of the outflow (e.g., the development of 
shocks or the initiation of Fermi acceleration) , it is possi- 
ble that high energy radiation such as gamma-rays may 
be produced in ways that make it possible to discriminate 
between single and binary black holes, or even possibly 
between slowly-rotating and rapidly-rotating black holes. 

In summary, our exploratory calculations represent a 
step towards realistic simulations of the accretion flow 
near binary supermassive black holes close to merger. 
Under certain conditions our results show differences be- 
tween single and binary black holes, and between non- 
rotating and rapidly rotating binary black holes, which 
is encouraging for future observations. As such, they 
motivate more realistic hydrodynamic and magnetohy- 
drodynamic calculations, which will be our next steps. 
If some of the qualitative differences found in our test 
particle calculations (with isotropically distributed ve- 
locities) are maintained in future work with fluids, this 
will imply prompt electromagnetic signatures that can be 
correlated with gravitational radiation signals detected 
by the space-based Laser Interferometer Space Antenna 
(jBaker et al.ll2007t iLang fc Hughes! I2009D . allowing pre- 
cise cosmolqgical probes and constraints on dark energy 
(|Arun et al.ll2009[ ) and even testing fundamental princi- 
ples such as the relative speed of propagation of photons 
and gravitational waves. 

We performed these calculations on Discover at NASA 
/ GSFC and Pleiades at NASA / AMES. GSR and MCM 
acknowledge partial support from the National Science 
Foundation under grant AST 06-07428. MCM was also 
supported in part by NASA ATP grant NNX08AH29G. 
The work at Goddard supported in part by NASA grant 
06-BEFSG6-19. JHW, STM and BJK were supported 
by appointments to the NASA Postdoctoral Program at 
the Goddard Space Flight Center, administered by Oak 
Ridge Associated Universities through a contract with 
NASA. 
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